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ABSTRACT 



A methodology is developed which determines the optimal 
allocation of patrol forces among selected deployment sites. The 
procedure uses a linear programming algorithm which minimizes a 
linear cost function, subject to restraining equations representing 
the total hours available, the relationship between on -station and 
transit hours, and base loading. A computer program is presented 
which translates input data into the format required by the IBM 
Mathematical Programming System/360 for the problem solution. The 
methodology can be utilized to determine the allocation of forces 
among selected bases, reallocation of forces when a base or bases 
must be removed from consideration, and the effect of utilizing 



additional bases . 



1 « 



DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 

j.val Postgraduate School MONTEREY, CA 93943-5101 

t-c.iterey, California 9394Q 

TABLE OF CONTENTS 



I. INTRODUCTION 11 

A. BACKGROUND 11 

B. OBJECTIVE 13 

C. ORGANIZATION 14 

II. SYSTEM DESCRIPTION 16 

A. AIRCRAFT 16 

B. BASES 17 

III. METHODOLOGY 2 0 

A. COSTS 22 

B. FORMULATION 2 6 

1. Notation 2 6 

2. On-Station Hours 2 7 

3. Transit Hours 2 7 

4. Total Hours Available 2 9 

5. Base Loading 29 

C. STATEMENT OF PROBLEM 3 0 

D. SOLUTION PROCEDURE 3 0 

IV. DISCUSSION 32 

A. EXTENSIONS 32 

B. ASSUMPTIONS 33 

C. RECOMMENDATIONS FOR IMPROVEMENT AND 35 

FURTHER STUDY 



i 



.. i v r* ' 

■ •• j v * 

/ ■> * 

V. SUMMARY 

APPENDIX A - RELATIONSHIP BETWEEN OPERATING RADIUS AND 

AVAILABLE ON -STATION TIME 

APPENDIX B - MAXIMUM FLIGHT-HOUR CAPABILITY 

APPENDIX C - UTILIZATION 

COMPUTER OUTPUT 

COMPUTER PROGRAM 

BIBLIOGRAPHY 

INITIAL DISTRIBUTION LIST 

FORM DD 1473 



40 

44 

47 

59 

68 

71 

72 

73 



4 



LIST OF TABLES 



Table 

I. 

II. 



Page 

Aerodromes of the Western Pacific Capable of 18 

Supporting P-3 Aircraft 

Annual Operating Costs 24 



5 



LIST OF ILLUSTRATIONS 



Figure 




Page 


1 . 


Possible Coverage of the Western Pacific by 
Patrol Aircraft Based in Selected Areas 


19 


2. 


Area A and Grid Overlay 


21 


3. 


Gross Weight versus Distance to On-Station 


42 


4. 


On-Station Hours versus Distance to On-Station 


43 


5. 


Pictorial Representation of Sample Problem 


48 


6. 


Sample Input Data for MPS/360 


51 


7. 


Control Language Program 


53 


8. 


Sample Computer Output 


57 



7 



LIST OF SYMBOLS 



(i,J) 


th 

= subarea formed by intersection of i — row 
and j — column 


XL 


= length of side of subarea parallel to x-axis 


YL 


= length of side of subarea parallel to y-axis 


(x k' Y k } 


= location of base k_ 


R... 

ijk 


= distance from base k_ to area (i,j) 


T 

OS k 


= available on-station time from base k 


T k 


= average sortie length from base k. 


CFH 

k 


= cost per flight -hour when flown from base k. 


C ijk 


= cost per on-station hour in area (i,j) when flown 
from base k_ 


b 

ij 


= on-station time required in area (i,j) 


X... 

ijk 


= on-station hours allocated to area (i,j) from base k. 


X m+i, jk 


= transit hours flown between area (i,j) and base k_ 


x . 
mis 


= total hours available for training and miscellaneous 
flying activity 


a k 


= flight time available from base k_ 


A 


= total flight-hours available 



9 



I . INTRODUCTION 



A. BACKGROUND 

At the present time the deployment concepts associated with the 
Navy's patrol aircraft in the Pacific Theater are little removed from 
those which evolved following the close of World War II. A majority 
of the advance bases currently supporting U. S. Naval Forces in the 
Western Pacific were acquired during the years following the Second 
World War. At that time the predominate thought concerning the 
positioning of advance forces was that the first line of defense should 
be as far away from the continental United States as possible. Covering 
nearly all of the transit routes between the Asian mainland and the 
Central Pacific, this chain of bases has provided the United States 
with a convenient surveillance platform. 

As long as the continued presence of the United States is re- 
quired in the Western Pacific to protect U. S. interests, the Navy must 
be ready to provide adequate forces for the following: 

1. Control of the sea-lanes and sea-areas against threats to United 
States interests, forces or commitments. 

2 . Continuing peacetime deployments in order to deter aggression and 
to support United States policy as it may evolve. 

3. Special surveillance, intelligence, and counter-surveillance 
operations . 
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It may be assumed that due to U. S. commitments established 



under the United Nations Charter, participation in SEATO and the ANZUS 
agreement, and many bi-lateral agreements and assurances that the 
advanced deployment of U . So Naval Forces in the Western Pacific will 
be required into the 1970's . 

Since naval forces are to be deployed during the next several 
years in approximately the same areas where they have been deployed 
over the past 10 years, the existing base structure may be regarded as 
adequate. It would be difficult to improve the geographical positioning 
of the present base structure without moving onto the Asian mainland, 
which is an alternative many military planners do not wish to consider. 

While the commitment of U . S. forces overseas is very likely to 
continue at or near its present level for the next few years, the con- 
tinued use of all present bases for the same time span is hr considerable 
doubt. It is entirely possible that continuing political pressure by 
groups in host countries may result in the denial of some bases to U. S. 
forces; for 'example, the Status of Forces Agreement with Japan is up 
for optional termination after 1970 on twelve months notice. 

Thought has already been given to a retrenchment to Guam, the 
only base site in the Western Pacific to which the U. S. has continuing 
access, and to the Micronesian Islands, which the U. S„ holds under 
a United Nations trusteeship. Called a "strategic trusteeship," it 
allows the U. S. to erect fortifications and garrison troops on the islands. 
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B. OBJECTIVE 



The increasing possibility of base denial and the rising cost of 
operating and equipping overseas forces have brought about the need for 
a reappraisal of present deployment concepts and the development of a 
method for the optimal allocation of available forces among available 
bases . 

It is the purpose of this thesis to present a method with which 
operational commanders may optimally allocate the patrol forces at 
their disposal, subject to operational requirements, operating areas, 
and forces available. 

The procedure developed requires as input data, information 
concerning the location of existing bases, the desired coverage of 
surveillance areas, and the amount of flight time available. Utilizing 
the Mathematical Programming System/360 Linear Programming package 
(MPS/360 LP) , available for the IBM 360 computers, a solution is 
determined which provides a minimum cost allocation of flight-hours 
among participating bases. 

The number of aircraft required at each location may be determined 
by comparing the number of flight -hours required with the flying hour 
capability of the aircraft. Since it is unlikely that this comparison 
will result in an integer solution for the number of aircraft required, it 
is necessary to round off to the next higher integer value. This will 
generally result in additional flight-hours being made available for 
training flights and other uses. Appendix B, combining the methods of 
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Sunde [1] and Mooz [2], presents a formulation for determining the 
flying hour capability of an aircraft from a knowledge of its operating 
hours and available maintenance data. 

C . ORGANIZATION 

In the formulation of this methodology the basic system considered 
is the P-3 series land based patrol aircraft and its supporting bases. 

No distinction is made between the various models of the basic P-3 
aircraft. 

A brief description of this aircraft, its operating characteristics, 
capabilities, and requirements is contained in Section II of this thesis. 
Also contained in Section II is a listing of some of the overseas bases 
capable of supporting P-3 operations. 

Section III presents the development of the methodology, A general 
linear programming formulation is followed in which a linear objective 
function denoting cost is minimized subject to a series of constraining 
relationships . 

Section IV discusses possible extensions of the methodology, 
inadequacies of some of the assumptions, and areas in need of further 
study. The thesis concludes with Section V, which presents a summary 
of the development. 
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Three appendices, A, B, and C, provide supplementary information. 



Appendix A contains the development of a linear approximation of the 
relationship between operating radius and on-station time . 1 Appendix B 
presents a method of determining the maximum flight -hour capability 
of an aircraft from available operational and maintenance data. In 
Appendix C, a sample problem is solved to demonstrate the use of the 
methodology. Also presented is the computer program, written in 
FORTRAN IV, which converts the input data for a problem into the format 
required for input into the linear programming algorithm. 



On-station time is defined to be that time spent in a specific 
operating area and does not include time necessary to transit to and 
from the operating area . 
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II. SYSTEM DESCRIPTION 



The system referred to in the section heading is considered to 
mean the P-3 series aircraft and its supporting bases. Although some 
earlier P-2 series aircraft are still in use, the fleetwide transition to 
the P-3 is sufficiently well along that only the P-3 will be considered 
in this thesis . 

A. AIRCRAFT 

The P-3 is a four-engine, low-wing, all-weather aircraft 
designed for patrol operations and antisubmarine warfare. It is in the 
12 7,000-pound gross weight class and is powered by four turboprop 
engines. The aircraft is fully pressurized and is capable of operating 
at all altitudes from Sea Level up to 34,000 feet and at speeds of from 
150 to 400 knots. As presented in Appendix A, during a normal mission 
time of 11.2 to 12.0 hours, the P-3 can transit to an operating area 
at a distance of over 1300 nautical miles and remain on-station for a 
period of four hours . 

The aircraft is normally manned by a crew of 12 men consisting 
of a pilot, copilot, navigator, tactical coordinator, flight engineer, 
and six technical specialists. 

Under normal operating conditions the aircraft will fly "profile" 
missions. Utilizing this "profile" concept, the aircraft will transit to 
a patrol area at altitudes between 17,000 feet and 22,000 feet at a 
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speed of 300-330 knots. The enroute altitude will generally depend 

upon the wind at different altitudes, distance to operating area, and 

takeoff weight. Upon arrival in the operating area, the aircraft descends 

to search altitude and reduces to maximum endurance airspeed. It is 

during this on-station period that one or possibly two of the aircraft's 

2 

engines may be "feathered" to increase the available on-station time. 
The return trip is usually made at a altitude of 2 5,000 feet to 30,000 
feet . 



B. BASES 

By considering the operating requirements of the P-3, the takeoff 
and landing distances, the fuel required, the necessary personnel, and 
the aircraft support requirements — and by referring to a listing of the 
major aerodromes is the Western Pacific, it is possible to compile a 
list of feasible operating bases for the P-3 aircraft. Table I presents 
a listing of bases which might be selected. 

Utilizing Table I and the information on operating radius versus 
on-station time as presented in Appendix A, Figure 1 may be drawn. 
From Figure 1 it can be observed that the P-3 aircraft, operating from 



A feathered engine, in this case, refers to one which has been 
shut down by the pilot to conserve fuel but which may be started at 
a later time . 
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suitable bases, can provide at least four hours of on-station coverage 
over a majority of the ocean area of the Western Pacific . It should be 
noted that in many areas a significant amount of overlap is provided. 

It is the optimal coverage of these areas of overlay which the methodology 
seeks to provide. 



TABLE I 

AERODROMES OF THE WESTERN PACIFIC 
CAPABLE OF SUPPORTING P-3 I ARC RAFT 



Tapan 


Misawa AFB 
Tachikawa AFB 
NAS Atsugi 
MCAS Iwakuni 


Okinawa 


Kadena AFB 
NAS Naha 


Guam 


Anderson AFB 
NAS Agana 


Philippines 


Clark AFB 
NAS Cubi Point 
Naval Station 
Sangley Point 


South Vietnam 


Danang 
Cam Rahn Bay 
Tan Son Nhut 


Taiwan 


Tainan 
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FIGURE 1 

POSSIBLE COVERAGE OF THE WESTERN PACIFIC 
BY PATROL AIRCRAFT BASED IN SELECTED AREAS 
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Ill . METHODOLOGY 



In the development of the methodology necessary for the optimal 
allocation of available resources, it will be convenient to assume that 
an area, A, exists into which it is desired to allocate a specified 
amount of patrol effort. This desired allocation will be measured in 
hours and will be assumed to constitute only on-station time. Located 
around, and within, area A are bases from which the required patrol 
effort is to be initiated. 

To facilitate the development, a rectangular grid will be super- 
imposed upon area A and its supporting bases such that the north-south 

V 

axis of A is aligned with the vertical axis of the rectangular grid. This 
grid is to be of sufficient size that all of area A and its supporting bases 
are enclosed within the borders of the rectangle. A Cartesian coordinate 
system is then established with the northwest corner of A as the origin, 
the positive x-axis lying to the east of the origin and the positive 
y-axis lying to the south of the origin. Distances along the coordinate 
axes will be measured in nautical miles utilizing the same scale as 
area A. The rectangular grid will subdivide area A into a number of 
subareas of equal size. The total number of subareas is the product 
of the number of columns (n) and the number of rows (m) within the 
rectangular grid. Assignment of a number _i, ranging from one to m to 
each row, beginning with the uppermost, and a number j_, ranging 
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from one to n to each column, beginning with the left hand side of A 
allows each subarea to be denoted by a pair of numbers, (i,j). Figure 2 
summarizes the development to this point. 




It is now possible to locate any point within the area enclosed by 

the rectangular grid by either of two methods. For example, the location 

of the point k in Figure 2 may be expressed as (2,5), indicating that it 

is within that subarea formed by the intersection of row 2 and column 5; 

or as (x, ,y, ) , which indicates that k lies x, miles to the right of and 
k k — k 

y^ miles below the origin. By choosing the spacing of the grid lines to 
be equidistant it is possible to assign a name XL to the length of the 
side of a subarea parallel to the x-axis and a name YL to the length 
of the side of a subarea parallel to the y-axis. 

For the purposes of this thesis it will be appropriate to assume 
that any flight designated to operate in a specific area will proceed 
to the center of that area prior to beginning its on-station period. 
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The distance, denoted R,,. , between any point k and the center 

ijk - 

of any specific subarea (i,j) may be written as a function of the 



coordinates of the point k and the location of subarea (i,j) in the 



following form: 



R 



[TiYL-^-y k ) 2 + (jXL - 



XL 



X 




2 I 1/2 



ijk 



2 



k 



If (x, ,y ) is in fact the location of base k, then R.., represents 

K. K 1 J K 

the distance in nautical miles from base k to operating area (i,j). 

As developed in Appendix A, the available on-station time from 



linear function of the distance between the base and the operating area, 
and the average sortie length, in hours, T^. 



Further utilizing the results of Appendix A, the maximum desirable 
operating radius, that which yields an on-station period of at least 
four hours, is found to be approximately 1350 nautical miles. 

A. COSTS 

In any problem requiring an optimal allocation of scarce resources 
it is necessary to evaluate the desirability of each possible alternative. 
By assigning a weighting factor, measured in dollars, to each variable, 
it becomes possible to express, in consistent terms, the value associated 
with each relationship. In the allocation of flight-hours, and hence 



base k, T /in any subarea, per sortie, may be approximated by a 
os, 



T 



os 



k 



T, - 0.0052R... 
k i]k 
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aircraft, among available sites it is desirable that this factor reflect 
differences in operating conditions, geographical relationships, and 
the level of operations . 

The system under consideration, that of patrol aircraft and bases, 
has been in the operating forces for many years. It is not required to 
consider any costs which might have been associated with any Research 
and Development, or Investment phase. The annual operating costs, 
those recurring outlays which are needed to operate and maintain 
activities in service, the only costs which need to be considered. 

Large [3] presents the listing shown in Table II, representing a partial 
breakdown of annual operating expenses. 

Examination of those areas listed in Table II discloses several 
which may be omitted from consideration. PAY AND ALLOWANCES are 
not directly related to the number of flight hours. Service personnel will 
be paid whether or not they fly. Similarly, TRAINING and ADMINISTRATIVE 
AND SUPPORT COSTS must be met even when no flying is performed. 

Items which do lend themselves to this type of consideration as a direct 
reflection of flying activity include, FUELS, LUBRICANTS, AND CON- 
SUMABLES as well as some of the MAINTENANCE categories. Consumable 
items whose usage rates are directly attributable to flying activity 
include flight clothing, and expendable stores such as sonobuoys , 
underwater sound signals, and smoke lights. The repair rate for many 
"Black Box" items is closely related to flight activity. Unfortunately, 
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TABLE II 



ANNUAL OPERATING COSTS 

I. EQUIPMENT AND INSTALLATIONS REPLACEMENT 

A. Primary Mission Equipment 

B. Specialized Equipment 
C . Other Equipment 

D. Installations 

II. MAINTENANCE 

A. Primary Mission Equipment 

B. Specialized Equipment 

C . Other Equipment 

D. Installations 

III. TRAINING 

IV. PAY AND ALLOWANCES 

V. FUELS, LUBRICANTS, AND OTHER CONSUMABLES 

VI. SERVICES AND MISCELLANEOUS 

A . Tran s porta tion 

B. Travel 

C. Miscellaneous 

VII. ADMINISTRATIVE AND SUPPORT COSTS 
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the Navy does not have a satisfactory method of assigning a cost to 
the repair of a particular radio, radar, or other "Black Box" component. 

It therefore becomes impractical to include repair costs of repairable 
components in a cost which relates to flying activity. 

By comparing the total cost of fuel, lubricants, and consumable 
items required to operate for a specified period of time with the number 
of flight-hours flown during the same period it is possible to determine 
an average cost per flight-hour, denoted CFH . Determining this figure 
for each location will provide a measure of the cost of operating as 
influenced by geographical location, operational requirements, and 
local operating practices. 

This figure will now be utilized to develop a costing procedure 
which can be used for the comparison of selected alternatives. If CFH^ 
is the cost per flight-hour when flown from base Jc, then the cost of 
one hour of on-station time in any subarea (i,j) that may be reached 
from base k can be determined. 



COST PER ON-STATION HOUR = 



TOTAL COST OF FLIGHT 

MAX. NO. OF ON -STATION HOURS 



which yields , 



C <T k > (CFiy 

1J T, - 0.0052R... 

k 13k 
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or 



C 



ijk 



CFH, 

k 

0.0052R , 
ilk 



C... denotes the cost per on-station hour in subarea (i,j) when flown 
ijk 

from base k. 



B. FORMULATION 



Under the assumption of a cost function which has a linear 
relationship with the on-station hours, the flight-hour allocation 
problem may be formulated as one which may be solved with the pro- 
cedures of linear programming. The problem becomes one for which it 
is desired to fulfill the operational requirements in each subarea at a 
minimum cost subject to certain restraining conditions expressible 
as linear equations. 

1 . Notation 

Prior to a formal statement of the problem, notation must 
be established. If (i,j) denotes a particular operating area and k a 
specific base , where i=l,...,m, j = 1 , . . . , n , and k = 1 , . . . , p , 
then the following definitions will apply: 

ijk number of on-station hours per month allocated to 

area (i,j) from base k_ 

x . number of transit hours per month to area (i,j) from 

m=i , jk 

base k in support of x. 

ijk 
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C cost per on-station hour in area (i,j) when flown 

ijk 

from base k 



b.. on-station hours per month required in area (i,j) 

i] 

x , total hours available per month for training and 

mis 

miscellaneous flying at all bases 



a, flight time in hours per month available from base k 

k ~ 



In the flight-hour allocation problem it is necessary to allocate 

an amount, x.., , of on-station hours per month from each of p bases 
ijk 

among mn operating areas where C... is the cost of one hour of on- 
ijk 

station time in area (i,j) when flown from base 1c. Each operating area 

requires b. . hours of on-station time per month, 
i] 

The objective function, which represents the cost of providing 
the required on -station hours, may be expressed as. 



C = 



n m t 

^ f * C... x... . 
j = l i=l ijk ijk 



C is now to be minimized subject to the constraints presented below. 

2 . On-Station Hours 

On-station hours allocated to each area from all bases will 
equal the on-station hours required in each area. This may be written 
as 





for i = 1 , 
and j = 1 



m 

,n. 



3 . Transit Hours 

In the determination of the total number of flight-hours to 
be allocated from each base it is desirable to know the number of transit 



27 



hours necessary to provide the required number of on-station hours. 

Where R ^ is the distance from base k to area (i, j) the relationship 

between the on-station time and the transit time may be obtained. 

From Appendix A, the tradeoff between the on-station time and the 

transit time on an individual sortie has been shown to be 

T = T, - 0 . 0052R. ,, . 
os k ljk 

K, 



If x , is the number of on-station hours allocated to area 
ijk 

(i,j) from base k and T is the average on-station time per sortie, 

° S k 

the number of sorties flown may be described as 

x. 



NUMBER OF SORTIES = 



ilk. 



os, 



Similarly, if x , ,, is the average number of hours of 
m+i , jk 

transit time allocated to area (i,j) from base k, the average transit 
time is 

T = 0 . 0052R. . 

tr k 1,k 

The number of sorties flown is then, 



NUMBER OF SORTIES = 



X m+i , jk 
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Equating these two equations, the number of on-station 



hours may be expressed as a function of the number of hours spent in 
transit. 

x. x 

ilk = m+i , ik 

T 0 . 0052R. 

os ljk 

K. 

As determined previously 

T = T - 0.0052R. „ 
os k 1 jk 

which is substituted into the equation directly above, yielding, as a 
constraint; 



X i jk X m+i,jk 



- 1 



0.0052R 



ijk 



= 0 , 



4 . Total Hours Available 



The sum of all flight-hours allocated, including training, 



must equal the total hours available. 



p n 2m 

Si ,St A x ‘ jk Xmls 



= A 



The upper limit of "2m" in the summation over i^ indicates that both 
the on-station hours (i=l , . . . ,m) and the transit hours (i= m+1 , . . . , 2m) 
are to be added . 

5 . Base Loading 

The number of all flight-hours available at each base per 
month may or may not be known. If the capacity of a base is a 
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significant factor then an upper bound on the number of flight-hours 
available from base k may exist. If there exists an upper limit to the 
total available hours at any base k, this restraint may be expressed as: 



n 2 m 

2 . 2 x ^ a 

j=*L i=l ijk 1 



C. STATEMENT OF PROBLEM 

A complete analytical statement of the flight-hour allocation 
problem is now possible, to bring together the development of the 
preceding paragraphs. The problem is then to: 

p n m 

Minimize ,2, 2 .2, C... X... 

k=l 3=1 i=l ijk ijk 

for i=l,...,m, j = 1 , . . . , n , k = 1 , . . . , p 

Subject to , 



£ x. = b. . 
2 -., 13 k 13 



k=l 



x. - x 
13 k m+i, 3 k 



- 1 



0.0052R 



p A 

k=l j=l i=l 



ijk 

•” X 

ijk mis 



= 0 , 



Xj + x , _ = A 



n Zm ^ 

il. Ji Xi)k - 



D. SOLUTION PROCEDURE 

The linear programming problem formulated above is solved by 
the MPS/3 60 LP package through the use of a two-phase program in 
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which a routine written in FORTRAN IV translates the necessary input 
data into a format compatible with the MPS/360 LP requirements. When 
the transfer of input data has been completed, execution of the MPS/360 
LP portion of the program begins. A sample problem is presented in 
Appendix C and includes a discussion of the output from the MPS/360 LP. 
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IV. DISCUSSION 



A. EXTENSIONS 

Other areas to which the methodology presents an immediate 
solution concern the problem of base denial, the selection of alternate 
bases, and the problem of an increase in requirements after force levels 
have been established. 

The problem of base denial and the subsequent reallocation of 
forces may be simulated by removing a base from consideration in the 
problem formulation. This is readily accomplished by changing the ND 
entry on the data card for the appropriate base, as shown in Appendix C. 

The previously mentioned possibility of base denial raises the 
question of what alternatives are available if a base is lost. One 
solution is to reallocate available forces among the remaining bases with 
the hope of obtaining a feasible solution. Another is to consider the 
utilization of existing bases not presently supporting patrol forces, or 
the establishment of new bases. 

In any alternative which includes the introduction of a new base 
or the improvement of existing facilities, care must be taken to ensure 
that a detailed analysis of all requirements is made. It may evolve that 
it is less expensive to construct an entire new base than to provide 
for the incremental adjustments necessary to bring an existing base up 
to the capability required. Large [3] and WORC [5] have listed many of 
the items which must be taken into consideration. 
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One of the primary considerations in any comparison of alternatives 



is the effectiveness with which the requirements may be met. The 
methodology presented in this paper may be utilized to assist in this 
determination. By assigning an expected cost per flight-hour to each 
location, the alternate bases may be included in the flight-hour allo- 
cation procedure. In this manner the effect of each of the alternate 
sites may be observed. Objective results from the simulation may then 
be combined with the results of additional comparisons, both subjective 
and objective, prior to making the final decision. 

Requirements for a positive level of training hours or other flight 
activity may also be included in the solution procedure. If the require- 
ment is one covering all bases, the constraint, x , = b may be 

mis mis, 

placed into the program. To provide for separate requirements at 
selected bases, the constraint shown above must be broken down for 
each location, i.e., 

X misl ^misl' X mis2 ^mis2 * 

B. ASSUMPTIONS 

The formulation of the problem assumes that the total number of 
flight -hours available will be greater than the total requirement for on- 
station and transit time. If, however, the situation arises in which 
the requirements exceed the number of available flight -hours , additional 
procedures must be instituted. From an academic standpoint the problem 
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may be solved by the establishment of a fictitious base, a , , whose 

P+1 

available flight -hours are defined as the difference between the hours 
required and the total hours available., 

Written as a constraint this becomes: 

2m n 
2 
j=i 

The costs associated with the on-station hours flown between this 
fictitious base and each operating area should be related to the cost 
of being unable to furnish the desired coverage of the area. If such 
a quantitative figure cannot be determined, a cost of zero may be 
assumed which will then allocate flight -hours on a minimum cost basis 
to as many areas as possible. In actual practice the problem may be 
overcome by first comparing the total flight-hours available with those 
which result from an infeasible solution to the linear programming 
problem, A subjective decision must then be made as to the necessity 
of coverage in each subarea, and the amount of coverage desired. By 
reducing the total requirements a feasible solution to the problem may 
be obtained . 

The manner in which non-feasible base-area combinations are 
removed from consideration is in need of revision. A more positive 
method, rather than the assignment of high costs, is necessary. It is 
possible, in some circumstances, for an undesirable base-area 




x - x 
k=l ijk ijp+1 



= A 



n 

£ 

j=i 



2m 

£ 



x. - A 
ijk 



P+1 



xi 1 

k=l 
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allocation to enter the solution. Such a condition might arise during 
the solution in the case where the base nearest the area concerned is 
at its upper bound, if one exists, and all remaining bases are outside 
the operating radius of the aircraft. In this case, the solution pro- 
cedure will utilize the only cost available, $999, to achieve a minimum 
cost allocation . 

C. RECOMMENDATIONS FOR IMPROVEMENT AND FURTHER STUDY 
The procedure suffers from its dependence upon estimates of 
operational requirements. While it is possible to obtain objective 
values based upon past requirements, care must be taken to ensure 
that the figures are not inflated by subjective estimates of future 
requirements. An overestimation of these requirements, while providing 
an excess of available flight -hours for training purposes and unexpected 
demands, will result in a lower utilization of aircraft and flight crews. 
The rapid response capability of the P-3 (it is possible to position an 
aircraft and crew at any point in the Pacific within 24 hours) indicates 
that operational commanders should position their patrol forces at 
overseas bases such that the expected level of requirements is met. 
Unusually heavy and unexpected demands upon the system may be 
handled by releasing forces from their home port. An alternative method 
might be a probabilistic interpretation of the flight-hour requirements. 
This would enable the requirements to be structured such that any 
chosen level of operations might be handled. 
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The problem as stated does not take into consideration the 



possibility of a minimum acceptable level of operation at each base. 

If a minimum level does exist it may be inserted into the program by 
selection of an appropriate a, value and utilization of a greater-than- 

K. 

or-equal-to constraint relationship. 

An area which requires considerable study is that of the role 
played by the training requirements of a deployed squadron. Under the 
present structure, patrol squadrons are in a state of continual change, 
with deployed units being made up of both trained and partially trained 
personnel. This requires a continuing, heavy, training program which 
often suffers under the weight of operational requirements . Training 
needs on deployment are filled as the opportunities arise but are con- 
tinually outpaced by operational demands. It would appear that a more 
feasible approach to this problem would be the creation of a larger basic 
unit than the present squadron, which could then deploy a majority of 
trained personnel, reducing the training requirements at deployed sites 
to a minimum. 

Costs, though they continually play a large role in any problem 
related to the optimal allocation of resources, are among the more 
difficult items to identify. The expansion of the concept of a cost per 
flight -hour to include specific costs for operational, training, and the 
other types of flying performed, would greatly enhance the capability 
of the methodology by allowing a more complete breakdown of the 
requirements . 
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The assumption of a linear cost function should also be investi- 



gated. It is possible that the further division of the cost per flight- hour 
concept would result in the determination of a non-linear variation 
between the cost of operating in an area and the time spent in that area. 
Variables which might enter into the determination of a non-linear 
relationship include the type of search performed, weather, and search 
stores expended. 
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V. SUMMARY 



A method has been developed by which force commanders may 
optimally allocate the patrol forces at their disposal. This is accom- 
plished subject to operational requirements, operating areas, and the 
forces available. Provided input data defining the location of existing 
bases, desired coverage of surveillance areas, and available flight- 
hours, the methodology utilizes the Mathematical Programming 
System/360 to develop a minimum cost allocation of available forces. 
The number of aircraft required at each location may be determined by 
comparing the number of flight hours required within the flying hour 
capability of the aircraft. 

The inputs required for the computer formulation are, the on- 
station hours required in each subarea, the location of bases under 
consideration, the flight-hours available at each base, the average 
sortie length in hours, and the average cost per flight-hour for the 
aircraft. 

The outputs generated are, the total flight -hours required from 
each base, a complete breakdown of the on-station and transit hours 
flown from each base, the total time available for training and other 
missions, and the total cost of providing the on-station coverage 
required . 
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The methodology presented in this paper derives a large measure of 



its usefulness from its inherent flexibility. The sample problem, which 
consisted of 42 subareas and four bases, required a linear program with 
215 row constraints and 337 columns. The MPS/360 LP is capable of 
solving a linear programming problem with over 4000 row constraints 
and an unlimited number of columns. 

Alternate bases may be included in, or removed from, the solution 
procedure with a minimal amount of effort, thus providing a rapid, 
efficient, means of determining the role of each location in the overall 
picture . 

An increase requirement in any area after forces have been deployed 
may be handled by changing the required on-station time in the area 
concerned, and adjusting the a values of each base to reflect the 

K 

number of aircraft at each location. The methodology will then determine 
any necessary reallocation of forces to handle the additional requirements. 
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APPENDIX A 



RELATIONSHIP BETWEEN OPERATING RADIUS 
AND AVAILABLE ON -STATION TIME 

In determining the relationship between the operating radius and 
the on-station time per sortie it becomes convenient to make the 
following assumptions regarding the initial configuration of the aircraft: 

1 . P-3B , takeoff weight of 12 7, 500 pounds . 

2. Full fuel load of 59,800 pounds and 300 pounds of water, 

3 . Outbound flight at 18,000 feet to 22,000 feet altitude . 

4. Return flight at 28,000 feet, 

5. Zero-fuel weight of 67,400 pounds. 

6. Reserve fuel of 8500 pounds. 

7. Flight to and from the operating area will be flown according to the 
maximum range speed schedule as presented in the P-3A/P-3B Natops 
Handbook . 

Based upon the previously stated assumptions and utilizing the 
material in the P-3 Natops Handbook [9], Figures 3 and 4 can be con- 
structed. Figure 3 depicts the relationship between the gross weight 
of the aircraft, operating radius, and available on-station time. 

Figure 4 illustrates the linear relationship which exists between the 
operating radius and the available on-station time. 



40 



A least squares regression analysis of the sample points in 



Figure 4 results in the relationship 

T =11.2- 0.0052R 
os 

between the operating radius and the available on-station hours per 

sortie. T is the available on-station hours per sortie, R is the 
os 

operating radius in nautical miles, and 11.2 is the average sortie 
length in hours . 

Neither Figure 3 nor Figure 4 takes into account the increase in 
on-station time possible if one or two engines are feathered. The 
estimates may therefore be considered to be slightly conservative 
and more useful for planning purposes. 
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DISTANCE TO ON -STATION 
(NAUTICAL MILES) 



APPENDIX B 



MAXIMUM FLIGHT HOUR CAPABILITY 



The lifetime of an aircraft can be divided into a combination of 
flying time and ground time. Flying time can be broken down into 
separate categories to indicate the type of flying performed. Examples 
of these might be; (1) operational, (2) training, (3) repositioning. For 
the purpose of determining the maximum flight hour capability all flight 
time can be treated the same. 

Ground time can be divided into the following divisions; (1) Ready- 
alert and standby, (2) undergoing maintenance, (3) awaiting spares, 

(4) turn-around time, (5) operationally ready but not flying. In keeping 
with present Naval terminology (2) and (3) will be referred to as, 

(2) not operationally ready due to maintenance (NORM) and (3) not 
operationally ready due to supply (NORS) . 

The total number of hours available for flight per month per aircraft 

i 



(average) is 730 hours, as determined by: 

24 (hours/day) x 365 (days/year) 



HOURS PER MONTH = 



12 (months/year) 
HOURS PER MONTH = 730 hours/month. 
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These 730 hours of available time per month per aircraft may be 



grouped as follows: 



F 


FLIGHT HOURS 


GA 


GROUND ALERT HOURS* 


GM 


NORM HOURS 


GS 


NORS HOURS 


GT 


TURNAROUND HOURS 


GO 


OPERATIONALLY READY BUT 
UNSCHEDULED 


D 


OTHER 


* includes ready-alert and standby 



For further development of the maximum flying hour capability of 

the aircraft it will be necessary to determine the number of NORS and 

NORM hours per flight hour. The number of NORM hours per flight hour 

for each aircraft may be determined in the following manner. Let be 

the number of NORM hours per flight hour, then 
GM 



Similarly K , the number of NORS hours per flight hour is found to be, 
s 



GS 




The number of available flying hours per month can now be seen to 
be limited by that time which must be allocated to maintenance, awaiting 
spare parts and other ground activities. These limitations may be 
expressed analytically as follows, for each aircraft 

F + GA + GM + GS + GT + GO + D = 730 hours/month 

and since GM = K F and GS = K F 
m s 

F( 1 + K + K ) + GA + GT + GO + D = 730 hours/month 
m s 
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which yields 



730 - ( GA + GT + GO + D ) 

^max 1 + K + K 

m s 

By minimizing or eliminating the time an aircraft is "operationally 
ready but not flying" ,(GO) , and those unexplained hours, (D) , this 
equation will establish the maximum flying hour capability of the aircraft 
consistent with current maintenance practices. 
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APPENDIX C 



UTILIZATION 



Utilization of the previously developed methodology will now be 
demonstrated by applications to a sample problem. Following the for- 
mulation of the problem; the computer program, preparation of the 
required input data, and the information contained in the computer output 
will be presented. 



A. SAMPLE PROBLEM 

Assume that the operating area is positioned as illustrated in 
Figure 5. The grid overlay has subdivided the area into 42 subareas, 
six rows and seven columns. Each subarea is assumed to be 300 miles 



on a side, yielding a total area covered of 1800 miles by 2100 miles. 



The four bases shown have the following coordinates, relative to the 



origin of the grid: 

Base X-coordinate Y-coordinate 



1 

2 

3 

4 



1020 

1140 

2040 

350 



1060 

240 

480 

1380 



The arcs around each base indicate the maximum practical operating 



radius for that base. 



The mission of the patrol forces assigned to these bases will be 
to provide coastal surveillance coverage of specific areas as indicated 
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FIGURE 5 

PICTORIAL REPRESENTATION OF SAMPLE PROBLEM 
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by the straight line segments in the figure. Additional requirements 
dictate the need for additional coverage in adjacent areas. From a 
knowledge of the type of forces available and the operational require- 
ments it is possible to estimate the on-station hours required in each 
subarea for a specified period. Assume that this has been done for a 
period of one month and is indicated by the small numbers in each box. 

If a subarea contains no number indicating the requirement, a requirement 
does not exist. 

The total requirement for on-station hours in the sample problem 
is then 2 900 hours per month. Assume now that the total number of 
flying hours available in this area per month will be 5500 hours. This 
figure includes, on-station hours, transit time, training hours and any 
other flight time. 

Base 4 will be assumed to be in an overloaded status and capable 
of supporting only a limited number of aircraft for patrol purposes . This 
will be indicated by placing an upper bound on the number of flight -hours 
available at base 4 of 600 hours per month. The remaining bases, 1,2, 
and 3 are capable of handling any number of aircraft that might be 
expected . 

Appendix A indicates that the average sortie length utilizing the 
"profile" concept for maximum aircraft utilization will be approximately 
11.2 hours. For the sample problem, assume that this figure will apply 
to each base. 
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The last figure required is that of a cost per flight-hour, CFH. 



This cost may be expressed as its true value, or as a multiple of a 
base value. For example, if the cost per flight -hour figures for bases 
1 through 4 are: $28, $33, $36, $31, they might be also presented as 
multiples of one of the values, say $2 8. In this form they would be 
presented as 1.000, 1.178, 1.391, and 1.107. 

B. COMPUTER PROGRAM 

The computer program performs the function of translating system 
requirements into the form required by the MPS/360, then executing 
the linear program and obtaining an optimal solution to the problem. 

The program consists of two parts, a routine written in FORTRAN IV 
which formulates the input data required for the MPS/3 60 and places it 
into storage. An MPS/360 program which retrieves the input data from 
its storage location, initiates a linear programming solution procedure 
and determines the optimal allocation. 

Inputs to the FORTRAN program are discussed below. After 
receipt of the input data the FORTRAN program computes the cost per 
on-station hour utilizing the relationship developed in Appendix A. If 
the range to any area is found to be greater than the maximum desirable 
operating radius of 1350 nautical miles, a cost per on-station hour of 
$999 is assigned to forestall inclusion of a non-feasible base-area 
combination. If an operating area lies outside the range of all bases 
considered in a particular problem, the requirement for that area is re- 
duced to zero, removing it from consideration. 
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The routine then computes the data required by the MPS/3 60, 



placing it into storage in the proper sequence. Figure 6 presents an 



example of the type of data and format necessary for the input to the 



MPS/360 program. 

NAME FLTHRS 

ROWS 
N COST 
E R1 



E R2 5 
G R2 6 
L R2 7 
COLUMNS 



Xlll 


COST 


27.63 


Xlll 


R1 


1.00 


Xlll 


R16 


1.00 

• 


• 

• 

X353 


• 

• 

COST 


• 

13.26 


X353 


R1 


1.00 


X353 


R26 


-.63 


RHS 






B 


R1 


100.00 


B 


R2 


200.00 

» 


B 


R25 


6000.00 



ENDATA 

SAMPLE INPUT DATA FOR MPS/360 
FIGURE 6 
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The first card contains the data set name, FLTHRS, and the last 



card, ENDATA, signifies the end of the data set, ROW cards specify 
the name to be assigned to the rows of the linear programming matrix, 
as well as the type of constraint (equality, inequality, or no constraint) 
represented by the row, COLUMN cards specify the name to be 
assigned to the columns in the linear programming matrix, and define, 
in terms of column vectors , the actual values of the matrix elements . 

RHS cards are used to specify the name of the right-hand -side con- 
straint vector. They are also used to define, in terms of column vectors, 
the values of these elements. Referring to Figure 6, the following 
interpretations are made. In the ROWS section, "N COST" indicates 
that this is the row corresponding to the objective function of the 
problem and does not have a constraint. "E Rl" signifies that 

V 

row Rl is an equality constraint while for row R2 6 the constraint 
relationship is greater -than or equal-to. If the only elements in row Rl 
are found in columns XI 11, X353, and B, the first equation may be 
written as 

Xlll + X353 = 100.00. 

The remaining constraint equations to the problem are formulated in a 
similar manner. 

When the transfer of input data into storage has been completed, 
execution of the MPS/360 LP portion of the program begins. MPS/360 
is composed of a set of procedures, a subset of which deals only with 
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linear programming. The method of solution of the linear programming 
problem is the ordered execution of a series of these procedures. The 
user decides upon the method of solution and conveys this to the 
MPS/360 in the form of the MPS/360 control language. Figure 7 presents 
the control language program utilized for the solution of the flight -hour 
allocation problem. 



PROGRAM 

INITIALZ 

MOVE (XPBNAME , ' PBFILE ' ) 
MOVE (XDATA , ' F LTHRS ' ) 
MOVE(XOBJ, 'COST') 

MOVE (XRHS/B 1 ) 

CONVERT 

CRASH 

PRIMAL 

SOLUTION 

EXIT 

PEND 



CONTROL LANGUAGE PROGRAM 
FIGURE 7 

Complete information regarding the MPS/3 60 is available in Mathematical 
Programming System/360, (360-CQ-14X) Linear and Separable Programming 
Users Manual [4] . 

C . INPUT DATA 

Required data for the solution to the flight-hour problem is of 
three types: 
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1. Information regarding the size of the area involved. 

2. Flight-hour requirements for each subarea. 

3. Base locations, costs per flight -hour at each base, and 

base utilization. 

The data deck is made up of cards containing the above information in 
the order presented. 

1 . Area 

The first card of the data deck contains six numbers which 
relate to the number of rows and columns which make up the grid overlay, 
the number of bases in the area, the length of the sides of each subarea, 
and the total flight-hours available. This information is conveyed to 
the program by the following two cards which specify the order and the 
format of the data. 

READ (5,102) M, N,P, XL, YL, AVAIL 
102 FORMAT (315,3F10.0) 

For the sample problem, the input data for this section will appear 
as shown below", with the figures (x) indicating the column in which the 
first figure is placed. 

674 300. 300. 5500. 

(5) (10) (15) (18) (28) (38) 

2 . Flight-Hour Requirements 

The flight-hour requirements will be read into the program 
in an array of the same dimensions as the grid, utilizing the cards 
presented below. 
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READ (5,100) ( (B(I, J) ,J=1 ,N) ,1=1 , M) 

100 FORMAT (6F10.0) 

The sample problem will appear as follows, each line 
referring to a separate card. 







200. 


200. 


100. 


50. 








200. 


100. 


50. 


50. 






200. 


200. 


50, 


50. 






200. 


200. 


50. 


50. 


50. 








200. 


100. 


50. 


50. 






200. 


100. 


100. 


50. 


50. 






(1) 


(ID 


(21) 


(31) 


(41) 


(51) 



3, Base Information 

The last group of data cards specifies information about 
each base in the area. The cards; 

READ (5,101) (A (I) ,T (I) , CFH (I) ,X(I) , Y(I) , ND(I) , 

11=1, P) 

101 FORMAT (5F 10. 0,12) 

convey this information to the program. A(I) is a number which corresponds 
to the maximum number of flight-hours per month that a particular base 
is capable of supporting. If there is no expected limit this number 
will be zero. T(I) is the average sortie length, while CFH (I) is the 
cost per flight -hour. The cost per flight -hour may be represented in 
either of the two forms mentioned earlier but consistency must be 
maintained within the program. X(I) and Y(I) correspond to the location 
of each base. The last figure , ND(I) , represents base utilization and 
may be either zero or one. If a base is to be utilized in the solution 
procedure the number will be one, if the base is not to be utilized, 
zero will be used. 
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Returning to the sample problem, the last section of the data deck 
will consist of the cards shown below. 



0.0 


11.2 


28. 


1020. 


1060. 


1 


0.0 


11.2 


33 . 


1140. 


240. 


1 


0.0 


11.2 


39. 


2040. 


480. 


1 


600. 


11.2 


31 . 


3 50. 


1380. 


1 


(1) 


(ID 


(21) 


(31) 


(41) 


(52) 



D. OUTPUT INTERPRETATION 

Figure 8 represents a reproduction of several segments of the 
sample program output. The cost of supplying the required number of 
operational hours is found in the "ACTIVITY" column under the heading 
"SOLUTION (OPTIMAL)" to be $108757.55. The "(OPTIMAL)" indicates 
that an optimal solution was reached. Other possible results are 
" (NON -OPTIMAL)" and "(INFEASIBLE)." The next section, "SECTION 1 - 
ROWS," contains the activity levels of each row in the optimal solution. 
Rows Rl through R42 indicate the operational requirements in each 
subarea. Row R211 specifies the total number of hours available, 
while rows R212 through R215 indicate the total flight hours required 
at each base. The first line following the hours available corresponds 
to base 1, the second to base 2, and so on. The final section, 
"SECTION 2 - COLUMNS," provides a complete breakdown of the 
operational and transit hours flown between each area and each base. 

For example, column X321 indicates that base 1 is allocating 200.00 
hours of on-station time per month to area (3,2) and column XI 021 shows 
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SOLUTION (OPTIMAL) 

TIME = 3 .20 MINS. ITERATION NUMBER 23 1 

. . . NAME ACTIVITY . . . DEFINED AS 

FUNCTIONAL 108757.54745 COST 

RESTRAINTS B 



SECTION 1 - ROWS 



NUMBER 


. . .ROW. . . 


AT 


. . .ACTIVITY. . . 


1 


COST 


BS 


108757.54745 


2 


Rl 


EQ 


• 


3 


R2 


EQ 


• 


4 


R3 


EQ 


200.00 


5 


R4 


EQ 


200.00 


212 


R211 


EQ 


5500.00 


213 


R212 


BS 


2191.00 


214 


R213 


BS 


872.98 


215 


R214 


BS 


# 


216 


R215 


EQ 


600.00 


SECTION 2 


- COLUMNS 






NUMBER 


.COLUMN. 


AT 


. . .ACTIVITY. . . 


232 


X32 1 


BS 


200.00 


233 


X331 


BS 


200.00 


407 


X1021 


BS 


71.99 


408 


X1031 


BS 


7.17 


553 


XMIS 


BS 


1836.01 



SAMPLE COMPUTER OUTPUT 
FIGURE 8 
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that 71.99 hours of transit time are necessary to provide the 2 00.00 



hours of on-station time required in area (3,2). Since the index of 
transit time requirements runs from i = m+1 , . . . , 2m , column X102 1 
refers to the transit time to area (3,2) from base 1. The last entry in 
"SECTION 2 - COLUMNS," contains the total hours available for 
other activities, In the sample problem the value of XMIS is 1836.0 
hours . 
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COMPUTER OUTPUT 



SOLUTION (UPTIMAL) 

TIME = 1.25 MINS. ITfcRATICN NUMBER = 231 



. . . NAME. . . ... ACT I V ITY . . . 

FUNCTIONAL 108757. 54745 

RESTRAINTS 



OEFINED AS 

COST 

B 



59 



SECTION 


1 - ROWS 


NUMBER 


. . .ROW. 


l 


COST 


2 


R l 


3 


R2 


4 


R3 


5 


R4 


6 


R5 


7 


R6 


8 


R 7 


9 


R 8 


10 


R9 


11 


RIO 


12 


R 1 1 


13 


R 1 2 


14 


R 13 


15 


R 1 4 


16 


R l 5 


17 


R 16 


18 


R 1 7 


19 


R 1 8 


20 


R 1 9 


21 


R20 


22 


R21 


23 


R22 


24 


R23 


25 


R 24 


26 


R 2 5 


27 


R26 


28 


R27 


29 


R28 


30 


R29 


31 


R 30 


32 


R31 


33 


R 32 


34 


R 33 


35 


R34 


36 


R35 


37 


R36 


38 


R37 


39 


R 3 8 


40 


R39 


41 


R40 


42 


R41 


43 


R42 


44 


R43 


A 45 


R44 


A 46 


R45 


A 47 


R46 


48 


R47 


49 


R48 



AT ...ACTIVITY... 

eS 108757.54745 

EC 

EG 

EC 200.00000 

EQ 200.0000C 

EQ 100.00000 

EC 50.00000 

EC 
EC 

EC . 

EC 200.00000 

EQ 100.0C00C 

EC 50.C00CC 

EC 50.00000 

EG 
EC 

EC 200.00000 

EC 20C.COOCC 

EC 50.00000 

EG 50.CC00C 

EG 
EC 

EQ 200.000CC 

EC 200.00000 

EC 50.C000C 

EG 50.C00CC 

EC 50.00000 

EG 
EG 
EC 

EQ 2CC.OGOOC 

EQ 100.00000 

EC 50.00000 

EQ 50.00000 

EC 
EG 

EQ 200.C00CC 

EC 100.00000 

EQ 100.00000 

EG 50.00000 

EC 50.00000 

EG 
EG 
EG 



EG 

EG 

EC 
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NUMBER 


. . .ROW. . 


AT 


. . .ACTIVITY. . 


203 


R202 


EC 


• 


204 


R203 


EC 


• 


205 


R204 


EQ 


• 


2C6 


R205 


EC 


• 


2C7 


R 206 


EC 


• 


208 


R207 


EQ 


• 


2C9 


R208 


EC 


• 


210 


K209 


EQ 


• 


211 


R210 


EC 


• 


212 


R211 


EC 


5500. OOOCC 


2 13 


R212 


BS 


2191. 00C07 


214 


R213 


es 


872.98309 


215 


R2 14 


BS 


• 


216 


R215 


LL 


600.00000 
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SECTION 2 - COLUMNS 



NUMBER .COLUMN. 

217 X 1 1 1 

218 X 1 2 1 

219 X 1 3 1 

220 X141 

221 X 1 5 1 

222 X 16 1 

223 X 1 7 1 

224 X211 

225 X 22 1 

226 X 2 3 1 

227 X24 1 

228 X 2 5 1 

229 X261 

230 X271 

231 X 3 i 1 

232 X32 1 
2 3 3 X 3 3 1 

234 X 34 1 

235 X35 1 

236 X361 

237 X371 

238 X411 

239 X421 

240 X431 

241 X441 
292 X 4 5 1 

243 X461 

244 X471 

245 X 5 1 1 

246 X 5 2 1 

247 X 5 3 1 

248 X 54 1 

249 X 5 5 1 

250 X 56 1 

251 X 5 7 1 

252 X611 

253 X62 1 

254 X63 1 

255 X64 1 

256 X651 

257 X66 1 

258 X671 

259 X 1 1 2 

260 XI 22 

261 X 1 32 

262 X 142 

263 X 1 52 

264 X 162 

265 X 1 72 



AT 


. . . ACT IV ITY . . 


BS 


• 


LL 


• 


LL 


• 


LL 


• 


ES 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


200.C00CC 


ES 


• 


BS 


• 


LL 


• 


BS 


• 


BS 


• 


ES 


200.00000 


BS 


2CC.00CCC 


BS 


50.CCCCC 


ES 


50.C000C 


BS 


• 


LL 


• 


BS 


• 


BS 


200.000CC 


BS 


50.000CC 


BS 


5C . COOOC 


BS 


5O.CC0CC 


ES 


• 


BS 


• 


LL 


• 


BS 


196.23895 


BS 


ICC. COOOC 


BS 


50 . 000 OC 


BS 


5C. COOOC 


BS 


• 


ES 


• 


LL 


• 


BS 


• 


BS 


1CG.C0000 


BS 


50. CCOCC 


BS 


50.C00CC 


BS 


• 


BS 


• 


ES 


• 


8 S 


• 


BS 


200.0C00C 


BS 


20C.0C000 


BS 


ICO. coooo 


BS 


50.0000C 


BS 


• 
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NliMtiEK .COLUMN. 

266 X212 

267 X222 

2fcd X232 
269 X242 

27C X252 

271 X262 

272 X272 

273 X 3 1 2 

274 X 322 

275 X 3 32 

276 X 342 

277 X352 

2 7 d X362 
2 79 X 3 72 
280 X412 

26 1 X 4 2 2 

282 X432 

283 X442 

264 X452 

265 X462 

266 X472 

287 X 5 1 2 

288 X522 

269 X532 

290 X 542 

291 X552 

292 X 562 

293 X572 

294 X6 12 

295 X622 

296 X632 

297 X642 

298 X652 

299 X662 

3C0 X672 

301 X 1 1 3 
3C2 X 1 2 3 
3C3 X 1 3 3 
3G4 X 143 
3C5 X 1 5 3 
3C6 X 163 
3C7 X 1 73 
3C8 X213 

3C9 X223 

3 1C X233 

311 X 24 3 

312 X253 

313 X26 3 

314 X 2 73 

315 X 3 1 3 

316 X 3 2 3 



AT 


. ..ACTIVITY. . 


ES 


• 


es 


• 


8S 


• 


es 


ICO. 00000 


8 S 


50. COOOC 


ES 


50.00000 


BS 


• 


BS 


• 


8S 


• 


BS 


• 


BS 


• 


ES 


• 


BS 


• 


BS 


• 


es 


• 


ES 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 


BS 


• 


es 


• 


BS 


• 
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NUMBER 


.COLUMN. 


AT 


...ACTIV ITY . . 


317 


X3 33 


BS 


• 


318 


X 343 


BS 


• 


319 


X353 


BS 


• 


320 


X 363 


BS 


• 


321 


X 37 3 


BS 


• 


322 


X413 


BS 


• 


323 


X423 


BS 


• 


324 


X433 


BS 


• 


325 


X443 


es 


• 


326 


X453 


BS 


• 


327 


X463 


BS 


• 


328 


X473 


BS 


• 


329 


X 5 1 3 


BS 


• 


330 


X523 


es 


• 


331 


X 5 33 


BS 


• 


332 


X 54 3 


BS 


• 


333 


X553 


ES 


• 


334 


X 563 


BS 


• 


335 


X573 


BS 


• 


336 


X6 1 3 


LL 


• 


337 


X62 3 


BS 


• 


338 


X633 


es 


• 


339 


X643 


BS 


• 


340 


X 653 


BS 


• 


341 


X663 


es 


• 


342 


X 673 


BS 


• 


343 


X 1 1 4 


es 


• 


344 


X 1 24 


BS 


• 


345 


X 1 34 


BS 


• 


346 


X 144 


BS 


• 


347 


X l 54 


BS 


• 


348 


X 164 


BS 


• 


349 


X 1 74 


es 


• 


350 


X214 


BS 


• 


351 


X224 


BS 


• 


352 


X234 


8S 


• 


353 


X244 


BS 


• 


354 


X254 


es 


• 


355 


X264 


es 


• 


356 


X274 


BS 


• 


357 


X314 


es 


• 


358 


X324 


BS 


• 


359 


X334 


es 


• 


360 


X344 


BS 


• 


361 


X 3 54 


es 


• 


362 


X364 


es 


• 


363 


X374 


BS 


• 


364 


X414 


BS 


200.C000C 


365 


X424 


es 


• 


366 


X434 


BS 


• 


367 


X444 


BS 


• 



64 



NLNBER .COLUMN. 

368 X454 

369 X464 

370 X474 

371 X514 

372 X 524 

373 X 534 

374 X 544 

375 X 5 54 

376 X 564 

377 X 5 74 

378 X614 

379 X624 

38C X634 

381 X644 

382 X654 

383 X 664 

384 X674 

385 X 7 1 1 

3 86 X 72 1 
38 7 X 7 3 1 
383 X 74 1 

389 X 751 

390 X 76 1 

391 X 77 1 

392 X 8 1 1 

393 X 82 1 

394 X831 

395 X 84 1 

396 X 8 5 1 

397 X 86 1 

398 X871 

399 X911 

4 CO X92 1 

4C1 X93 1 

4 C 2 X941 

4C3 X95 1 

4C4 X961 

405 X 97 1 
4C6 X 10 1 1 
4C7 X 1 02 1 
4C8 X 10 3 1 
4 C 9 X1041 

410 X 1 0 5 1 

411 X 1 0 6 1 

412 X 1 0 7 1 

413 Xllll 

414 X 1 1 2 1 

415 X1131 

416 X 1 14 1 

417 X 1 1 5 l 

418 X 1 16 1 



AT 


. . .ACT IV IT Y . . 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


3.761C5 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


200.000CC 


BS 


100.00000 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


BS 


• 


LL 


• 


BS 


• 


BS 


• 


BS 


• 


LL 


• 


LL 


• 


LL 


• 


BS 


• 


BS 


• 


BS 


89.72633 


LL 


• 


LL 


• 


BS 


• 


LL 


• 


BS 


• 


BS 


86.2069C 


BS 


47. 18094 


BS 


8.45166 


BS 


13.30845 


BS 


• 


BS 


• 


LL 


• 


BS 


71.99424 


BS 


7. 17154 


BS 


.74503 


BS 


9.05141 


BS 


• 


BS 


• 


BS 


• 


BS 


82.87118 


BS 


22.54283 


BS 


7.82718 


BS 


12.81066 


BS 


• 
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NUMBER 


.COLUMN. 


AT 


. . . ACTIV ITY . . 


4 19 


X 1 1 7 1 


BS 


• 


420 


X 1 2 1 1 


BS 


• 


42 1 


X 1 22 1 


LL 


• 


422 


X 1 23 1 


BS 


43. 10345 


423 


X 1 24 1 


BS 


18.89645 


424 


X 12 51 


BS 


22.87283 


425 


X 1 26 1 


BS 


• 


426 


X 1 27 1 


BS 


• 


427 


X712 


BS 


• 


428 


X 722 


BS 


• 


429 


X732 


BS 


45.65168 


430 


X742 


BS 


12.56124 


431 


X752 


BS 


11.86662 


432 


X 762 


BS 


15.82779 


433 


X772 


LL 


• 


434 


X812 


LL 


• 


435 


X822 


LL 


• 


436 


X 832 


LL 


• 


437 


X842 


BS 


11.86662 


438 


X852 


BS 


7.99744 


439 


X 862 


BS 


17.2117C 


440 


X872 


LL 


• 


441 


X912 


LL 


• 


442 


X922 


LL 


• 


443 


X932 


LL 


• 


444 


X 942 


LL 


• 


445 


X952 


LL 


• 


446 


X 962 


LL 


• 


447 


X972 


LL 


• 


448 


X1012 


LL 


• 


449 


X 1 022 


LL 


• 


450 


X 1032 


LL 


• 


451 


X 1 042 


LL 


• 


452 


X1052 


LL 


• 


453 


X 1062 


LL 


• 


4 54 


X 1 072 


LL 


• 


455 


X 1 1 1 2 


LL 


• 


456 


XI 122 


LL 


• 


457 


X 1 1 32 


LL 


• 


458 


X 1 142 


LL 


• 


459 


X 1 1 52 


LL 


• 


460 


X 1 162 


LL 


• 


461 


X 1 1 72 


LL 


• 


462 


X 1 2 12 


LL 


• 


463 


X1222 


LL 


• 


46 4 


X 1232 


LL 


• 


465 


X 1 242 


LL 


• 


466 


X 12 52 


LL 


• 


467 


X 1262 


LL 


• 


468 


X 1 2 72 


LL 


• 


469 


X 7 1 3 


LL 


• 
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521 

522 

5 23 

524 

525 

526 

527 

528 

5 2 5 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

54C 

54 1 

542 

543 

544 

545 

54 6 

547 

548 

549 

550 

55 1 

552 

553 



.COLUMN. 


AT 


... ACT IV ITY . . 


X 844 


LL 


• 


X 854 


LL 


• 


X 864 


LL 


• 


X874 


LL 


• 


X914 


LL 


• 


X924 


LL 


• 


X 5 34 


LL 


• 


X944 


LL 


• 


X954 


LL 


• 


X964 


LL 


• 


X974 


LL 


• 


X 1 0 14 


6S 


43.64906 


X 1024 


LL 


• 


X1034 


LL 


• 


X 1 044 


LL 


• 


X 10 54 


LL 


• 


X 1064 


LL 


• 


X1074 


LL 


• 


X 1 1 1 4 


8 S 


• 


X 1 1 24 


BS 


. 1916C 


XI 134 


LL 


• 


X 1 144 


LL 


• 


X 1 1 54 


LL 


• 


X 1 1 64 


LL 


• 


X 1 1 74 


LL 


• 


X 1 2 1 4 


BS 


36.96858 


X 1 224 


BS 


15.42972 


X 1 234 


LL 


• 


X 1244 


LL 


• 


X 1 2 54 


LL 


• 


X 1 264 


LL 


• 


X 1 2 74 


LL 


• 


XMIS 


es 


1 836 . C 16 84 
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CC M PUT C R PROGRAM 



RROG ram COMPIJTFR I N'MJT HAT A RFOHIREP 
PQR thf SOI M T I c N nc THf fl T GMT HOUR 
A LL DC A T I CM PROBI FM AND orfcpmts t H c r)\jt 
IN A FORMAT C PMR AT I RLE WITH THF INPUT 
RFOIJI RF M FNTS FOR THE IB W SYSTEM 36) 
MATHENA TIC AL PRCGRAMMI NO, SYSTEM. 

VARTARLF NAMES AND PROGRAM IN d MT c 



w MDMppR CF ROWS CF P£OU I RF MFNT C ARRAY 

V MUMRFR r P COLUMN? CF RFOUIPFMFNT? ARRAY 

R MIIMRRR rF RASFS 

XL I. E NT,TH C F HORIZONTAL S IO c CF ? U R A R F A 

YL I F N'G th C= VFRTICAL S I OF re SUPARFA 

P( I,J) ...FLIGHT HHjo RFOUI RFMFNTS °FP AREA 

Xl^) X -COOPT I MATE cf Ri,Rc location 

V(M Y— COOROJ NA T E C c RASE LOCATION 

AVA I L. . . . T0t 6 L FLIGHT H0U° ? AV ML A RL c 

T ( K ) AVFRAHF SORTIE LENGTH FROM =AC.H RASE 

C.FHI K ) . . .r OST F r R FLIGHT HOUR FQOM EACH RASE 
NP f K ) . . . , P E NPTF c IF RASE IS TO RF UTILIZED 



I DC fodva T( AF10. C) 

1C 1 c OPMA T( «1R,0 ,1 2) 

1C 2 FOPVA T( 31 S,3F!^ .C) 

2C'C FQP v/\ T( *NAMF • ,T15 , 'FLTHRS • ,TRP , • A* / • ROWS • , Tp O , ' A ' / T2, 
1 'N • ,T5 , 'COST • , T P r , • A » ) 

2C1 forma t( t 2 , • = *,T5 » 'R » ,T3 ,T8u, ’A ' I 
2C2 c OR W A T ( *C CL U M NS • ,T»r, , • A •) 

2 C 3 FORMA T( T2 , 'R • f TR , »R • f 13 ,T90 » • A* I 

2C4 FQRM£T(T5 , 1 X ' *?I2»T1 C , *CO$T • , T2 *S , F 1 J . 3 , TPC ♦ • A • ) 

2CE form? T( Te , * x t ,3i 2 ,T1 5 , »R • , 13 ,T3 1 , • 1 .Cj • ,TPwi, • A» ) 

2C6 FOPMA T( T5 , • X • ,3T2,T15 , »R • , 13 ,T2R,FlO .3 ,tpo , • A» ) 

207 FORMAT! *RHS ', T«C » *A • ) 

2C? forma T( T E , 'R • , T 15 , *P • ,13 ,T25 ,F10.3 ,Tp:j , • A • ) 

2C° FOP ma T( ’EA'OA ta * ,TRO , 'A • ) 

210 FORMA T(T2,'L', t R,'R' ,I3,TSC,'&' ) 

211 FOPMATITE, 'XMI R • ,T1R , «R * ,1 3 , T 31 , *1.00' ,TF0, M* ) 

2 3 1C FOPMA T( T2, "= • , TK , « R • ,T 2 ,TflC , • A • ) 

20 PO FORMA T( T5 , 'R ' ,ti = , tp . , ,T25 ,Fl ).3 ,TR0 , • A • ) 

P I MENS I CN R< 20,20 ,X(1C I ,Y(10) ,R(2C , 2D , 1C ) , T { 10 I , 
1CFHI 10) ,C (2C,2C ,1C) ,A (10) ,ND(K ) 

I N T EGF R P , R x 

OAT* 7EPC ,RX/C. C ,0/ 

R C A F>( *5,102) M .M , P, XL, YL, AVAIL 
REAP («5,1C3) < < P< I , J) , J=1 ,NI ,1=1 ,M) 

P^AO <3,1011 (A ( n ,T( ! ) ,CFH(I ) ,X( I) ,Y( II ,NOM), 1 = 1, R ) 
T 1=M*N*( R + l I +1 
T P P W = I 1 + P 
M 1 = M + 1 
M 2=2*M 

no ic k=i,r 
O n ic I = 1 , M 
no 10 J = 1 » N 
Xl= J 
Y 1 - I 

R Y = A P $( Yl*YL-( YL/2.C) -Y( KI I 
RX = ABR(Xl*XL-tXL/2.0)-X(K) ) 

»< I ,J,K)=S0RT(PX**2.+RY**2. I 

CCMRUTF CORT RfR ON STATION HOUR 

C(I ,J»K)=CFH(K) /(!.-( 0.005 2*R(T ,J,K ) /t (KII I 
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IF 3 A NO F TO P P F R AT I NO A RFA TS ORFA T ER 
T HA N 13*0 NAUTICAL ”!LES ASSIGN * COST 
HP ♦CRR.nr prp HPI ) R 

no 1« I = 1 t o 

! 0 ?X=NP( I M-&X 

no 12 T = l, w 
12 ,1 = 1, N 

L = 0 

no n k=i,d 
R X = p ( I , J , x ) 

IF A N A 0 F A CANNOT RE P C »CHED FPflM ANY 
BA ?c R F PIJC F TTC RE O') I R ^ M ENT TQ 7PR0 

11 IF(o X .GT. 13SC. ) L = L+ 1 

12 IFfL.OF.oxj R( T , J) =0.0 
W° T TF ( o,20A) 

C P MD l TC Data fcp M p CW C " SEC T IOM 



DO l c I = 1 » ! °CW 
r c ( II — T I 11 , 17,17 



13 


K = ! - I 1 
















r m & t * 


) 


) 


14,14 


,16 




14 


T F r MO ( 


X 


) 


. F Q . 0 


) 


GO 


TO 17 


1 F 


WO T ( 


0 


♦ 


203 ) 


I 








on ▼r 


1 


c 










If 


WQ t rr ( 


Q 


♦ 


21 ?) 


I 








TO 


1 


c 










17 


I ^ ( I-l 


r\ 

U 


( 


) 170 


t 


171 , 


,171 


1 7 C 


wo ! ( 


3 


♦ 


2011) 




I 






pn to 


1 


c 










171 


WP! TP ( 


p 


♦ 


2 C 1 ) 


I 






l c 


CONTI Nl! 


c 











rpvouTF r A T A FHR "CPLIIMNS” SFCTIPN 

WR I TF ( o,202 ) 

DO 20 K=l, c 
DP 2C 1=1, ^ 
on 20 ,1 = 1 ,N 
T 2=1 1 + K 

! 2 3 5= ( T - 1 ) *N+J 

I 2 '■ 6 A =K*v*N 4-( I-l ) *N+.I 

wo I tf ( 9, 204) I , J,K ,C ( I , J,K) 

WOttf(o,2Co) I,J,K,!2C5 
Woy-Ff Q,2C5) I,J,K,I20FA 
WP I TF ( 0 , 20 S ) I ,J,K,T1 
20 WP T " r F ( o,2C5) I , J,K,T2 

DO 30 K=1,D 

D° 3( I= V 1, M 2 

DO 30 .1 = 1, N 
I2=! 1+K 

I 20 5 A =K*w*n+ ( I-M- l ) *M+ j 

I m = r- w 

OX=i.OC-( T(^)/(C. ,C0 C 2*R(I M »,I,K) ) ) 

WR ! T F ( o,206) I ,J,K,I205A,CX 
W^TTf F,205) I , J , K , ! 1 
3C W° I TF ( ° , 20 5 ) T , J,K,I2 
WR ! tf ( P , 2 1 1 ) II 

COMPUTE CAT A FOR »PH$" SECTION 

WP T TF ( c , 21 7 ) 
no 4C I=I,M 
DO 4.2 ,1 = 1 , N 
I20T=( I-l )*N + J 
I F ( I ? C 7— ICO) 400,^01,401 
4lC WP t TC( o,20O0) I 2C7 ,B ( I , .)) 
r;o t 4P 

4C1 WP!TE( p^C 3 ) 1207 , BU ,J) 



69 



4' CDNTINIJP 

DO 50 K = 1 t ° 

DO 5P 1 = 1, v 
D n sr, .1=1, N 
I 2 07/' *K*MN+(1 -1 ) *N+ J 
IP( I2C7/-13C) 5 C> C , ‘j 0 1 ,501 
e CD W D T T F ( 3 , 2C °0 1 !2G7A,7FRP 
0° TO 5C 

= 11 w° T TC ( c , 2 C 0 ) I2 r '7A,7F 0 C 
=0 CONTI NUP 

W’TTF f a , 2 C 8 » II, A VAIL 
HO 70 K=l,o 
I 2= 1 1 + *< 

70 WPTTff o,2C«> I 2 , A ( K ) 

W° t T F ( a, 209) 

cjpo 

FMO 



/ /GO .PTOPFOO 1 DP UNIT = SYSPA»DSN=F, ^PS , DCB= ( R EC F V= FP , 9U< c I 7 P= 3 
// L c ’Fri = ( '0) ,DISP=(NFW,PASS) ,SPAC P =(TRK»(30*?),°LSF) 

//GO.CY^tm op * 



INPUT DATA DECK 



A 


7 4 3CC. 


300. 


5 30C . 




c.o 


0.0 


2 CO. 


20-:'. 


100. 


50. 


- • 


c.o 


C.O 


2 C 0 . 


10 ; . 


50. 


c c. 


: 0 /** 
‘ • V 


C.O 


20 0. 


20; . 


50. 


57. 


0. 0 


r.c 


20 0. 


20 J. 


5 0. 


c <" . 


*c. 


0. C 


o. 3 


C .0 


200 


100, 


so. 


50. 


0 • 0 


r\ 

• v 


200 


100. 


ICC. 


p 0. 


30. 


r . u 


0 .D 


C .0 


11.2 


23. 


10 20. 


1060 . 


1 


C .0 


11.2 


33. 


1 1 40 . 


2 40. 


1 


C . 0 


11.2 


30. 


2040. 


4 30 . 


1 


= ce. 


11.2 


31. 


35 C * 


1330 . 


1 


/ / S 2 


FXFf LTNPRGG 










//^PSl 


. 3/ VP DP D SN = 


*. SI. Go. 


F T 08FCC 1 


, D I 3 p= ( DLD, 


PASS I 


/ / ’1 D SI 


.SYSIN DD * 











CONTPCL LANGUAGE PROGRAM FOP wo$/3f C 



PP CGPA M 
INITIAL 7 

MO VP ( XPRNAMF , • PPF I LF * ) 
M 0 VF ( XDA T A , 'FLTHPS • ) 

VP V c ( X03 J , *C CST 1 ) 

M 0 VF ( XPHS, *3 • ) 

CONVERT 

cr tup 

CP A PH 
PP TM» L 

solution 
p x i t 
RPNP 



//MP S2. SYSPP I NT DO $ D ACF=(CYl ,6) 

/ /MD S2 . $Y S I N PC DSN=*.$1.G0.FTG3F301 ,DISP = (OLP,DFLF T F) 
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